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Abstract 

We consider a simple inflation model with a complex scalar field coupled to gravity 
non-minimally. Both the modulus and the angular directions of the complex scalar are 
slowly rolling, leading to two-field inflation. The modulus direction becomes fiat due to 
the non-minimal coupling, and the angular direction becomes a pseudo-Goldstone boson 
from a small breaking of the global U{1) symmetry. We show that large non-Gaussianity 
can be produced during slow-roll inflation under a reasonable assumption on the initial 
condition of the angular direction. This scenario may be realized in particle physics models 
such as the Standard Model with two Higgs doublets. 
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1 Introduction 



Cosmic inflation pj has been the cornerstone beyond the standard big bang cosmology. It 
makes the present observable universe homogeneous and isotropie on large scales and gives 
rise to the primordial perturbations [2]. These primordial perturbations, when reentering the 
horizon after inflation, become the seed for the subsequent structure formation that leads to 
the temperature anisotropies in the cosmic microwave background and the inhomogeneous 
distribution of galaxies. They are very well constrained by observations to foUow Gaussian 
statistics and to give nearly scale invariant power spectrum [3]. Although there have been 
many proposals for inflation models consistent with observations [1], it seems hard to distinguish 
between different models from the current observations. 

As we are getting more precise cosmological data from the ongoing and future experiments, 
e.g. from the Planck satellite |5], we will be able to constrain the models of inflation beyond 
Gaussianity and power spectrum in next few years. Especially, the possibility of observing 
primordial non-Gaussianity from inflation has recently drawn a lot of attention from both 
cosmologists and particle physicists [6] . Since most slow-roU inflation models predict negligible 
non-Gaussianity [7], the detection or nuU-detection of non-Gaussianity will be an important 
discriminator for inflation models. 

The inflation models with non-minimal couplings [8] have been receiving a renewed interest 
from the recent claim that the Standard Model (SM) Higgs field can support inflation provided 
that it is non-minimally coupled to gravity [H] . The merit of the Higgs inflation consists in the 
minimality that the inflaton is the Higgs boson in the SM with a single additional parameter, 
the non- minimal coupling of the Higgs doublet. However, matching to the COBE normalization 
for the resulting density perturbation, one needs a large non-minimal coupling, given that the 
Higgs self-coupling appears large to be consistent with the Higgs mass bound at LEP (TU]. It 
has been shown that the large non-minimal coupling gives rise to a unitarity problem in the 
gauge boson scatterings of the SM [11] so the perturbative expansion breaks down at unitarity 
scale without adding extra degrees of freedom [121 [13]. Nonetheless, the predictions of the Higgs 
inflation may be maintained due to the fìeld-dependent unitarity cutoff [14j and/or depending 
on a ultraviolet completion [121 [I5]- Therefore, the Higgs inflation seems to remain as an 
interesting proposai . 

In this paper, we consider a simple two-field inflation model with a complex scalar field. We 
assume that the non-minimal coupling is not large such that there is no unitarity scale much 
below the Planck scale. The potential of the modulus of the complex scalar becomes fiat at 
large field value due to the non-minimal coupling. Moreover, a small breaking of the global 
U (1) symmetry gives rise to a fiat potential for the angular mode of the complex scalar too. On 
the other hand, the slow-roll parameters of the angular mode remain small during inflation such 
that the slow-roll conditions for the angular mode are always satisfìed. At horizon crossing, the 
vacuum expectation value of the modulus is chosen to obtain the required number of e-folds 
while that of the angular mode is set in such a way that it gives a negligible contribution to the 
slow-roll parameter e. Very interestingly, we show that a sizable contribution of the angular 
mode to e at the end of inflation leads to non-Gaussianity large enough to be detected. 

This paper is organized as follows. In Section |2l we first present the inflation model with 
a complex scalar field non-minimally coupled to gravity and show how two-field inflation is 
possible. Then, in Section [3] we analysis the perturbations during inflation using the 6N for- 
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malism and compute the non-linear parameter /nl as well as the power spectrum and 
the spectral index n^. Finally, in Section H] conclusions are drawn. 



2 Model 



We consider the general action for a complex scalar field with dimension-4 interactions in the 
Jordan frame, 



R-\d,<p\'-vM), 



(1) 

(2) 



Here we note that ^ and A respect the U{1) global symmetry while a, k and ( do not. For 
simplicity, we choose ali the dimensionless couplings to be real. Then, the unbounded below 
constraints give X + k + (>0, X-k + C>0 and ^/{X+~CY'-^ + A - 3C > 0. 



Writing the complex field in terms of two real fields h and x 



h 
7=2 



(3) 



the originai action ([T]) becomes 
Cj 1 



V-9J 



:R + 
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2x 
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Changing to the Einstein frame by a conformai transformation with 

n'^ = l+^h'^ + ah^ cos 



2x 
/ 



(5) 



we have 

2 
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1 + 

2x 
/ 



6^ 



4 + a cos ( — 
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2x 



1 + 



4 + a cos ( — 



4fì4 



sm 



2x 
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X + K COS ( ^ 1 + C cos ( 



f 



f 



(6) 



Thus, for a non-zero a, there is a kinetic mixing between the modulus and angular directions. 
For large values of the modulus, the potential becomes fiat along the modulus as in Higgs 
infiation. Meanwhile, the potential for the angular direction may be too steep if the symmetry 
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breaking terms proportional to a, k and ( are of order unity. To see this more clearly, we take 
the limit 1/ a/^. Then, ([6]) becomes 



V-9E 



4- 



1 + 



sm 



2x 
/ 



Qa 



sm 



1 A + ^cos(2x//) + Ccos(4x//) 
4 [^ + acos(2x//)]' 



(7) 



Then, in the potential, the first term proportional to A provides a vacuum energy for infiation 
while the pseudo-scalar x dependent terms can satisfy the slow-roll conditions for k, ^ A 
and a ^ ^. The small U{1) breaking terms are technically naturai in the sense that for equally 
small U{1) breaking parameters, the loop corrections to them are under control. In this case, x 
can be another infiaton component. The potential for x is of a very similar form to the infiaton 
potential for a codimension-two brane moving in a warped background in six dimensions [T7] . In 
the brane infiation case, infiation is assumed to end when the moving brane hits the background 
bran^. On the other hand, in our case, the modulus is the main infiaton component, which 
finishes infiation when the slow-roll condition is not satisfied any more. 

From now on, for simplicity, we consider the case with a = so that the kinetic mixing is 
absent. Also we work in the Einstein frame and drop the subscript E. The results remain valid 
for the more general case with non- zero a as far as |a| ^ this case, (E]) takes a much 

simpler form. 



-R- 



-9 



29? 



1 + 



id.xf 



A + K cos 



2x 
/ 



+ C cos 



4x 
/ 



We can find the canonical radiai field ip from 

d(p 1 



6^2/^,2 

dh V ^ ^ fì2 



as 



Vl + e^-isinh"^ + 6^)h\ - \/2tanh^^ 



1 + ^(1 + 60/^' 

Let US consider two different regime of h. First, for 0<h'^l/^,(f^h/Q and 

1 



(9) 
(10) 



A + K cos ( y I + C cos \ — 



Ili 



Meanwhile, for h ^ 1/^, we have a non-trivial relation 

dh TT^/^' 



from which we obtain 



(12) 



(13) 



^But this is not always the case: see e.g. Ref. [18] 
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Then, in this limit, the potential is approximated by 



4^ 



1 + (5i cos 



2x 
/ 



+ 82 cos 



4x 
/ 



(14) 



where 5i = k/\ and 82 = C,/\. Therefore, we have included a small potential terni for the 
modulus as well so that both the modulus and the angular directions of the complex scalar 
are slowly roUing. In order to canonically normalize the kinetic terni for the angular direction 
X during inflation, we need to choose / = 1/a/^, viz. the decay Constant of the would-be 
Goldstone boson is l/y^- This is comparable to a similar observation in Higgs inflation that 
the gauge boson mass during inflation is of (9(1/ y^), which is identified with the ultraviolet 
cutoff during inflation [13]. 

Focusing on the regime with /i ^ 1/^, we can write the Lagrangian as 



C 



(15) 



This is a particularly simple Lagrangian [19] which has been studied in models motivated from 
string theory: it includes a non-canonical kinetic term with 



and a potential of product form W{ip, x) = U{ip)V{x) with 



2^ 



(l-e-^WVey 



V{x) =1 + 6, cos (2 + à2 cos (4v^x) • 



(16) 

(17) 
(18) 



As will be shown later, the above product form of the potential makes it particularly easy 
to study two-fìeld inflation, not only for the background evolution but also for the explicit 
computations of the perturbations using the 6N formalism. 



2.1 Slow-roll parameters 

Since we are interested in the slow-roll phase of inflation, we can write the slow-roll parameters 
based on the potential derivatives. The first slow-roll parameter e can be written as 



1 ru'V e'^' rv 



2[u I ^ 2 [v 



(19) 



where a prime on U and V denotes a derivative with respect to ip and respectively. We can 
also define the second slow-roll parameter t] as 



V" 

V -e y 



(20) 
(21) 
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Likewise, let us define 



e'' =86'^ , (22) 
if =166" , (23) 

where a prime on h denotes a derivative with respect to y?. 

With the potential given by (ITTI) and (ITSl) . these slow-roll parameters are given by 

e''=o 2' (24) 

3 (i_e-2^/V6)2 ' 

2e [^isin (2v/ex)+2^2sin^ (4Vex)]' ^25) 

1 _ e-2^/v/6 [1 ^ (2v/ex) + ^2 cos (4v/ex)] ' ' 

= _ %-2^I^^^J± (26) 

3 (i_e-2WV6)2' ^ ' 

4^ ,5i cos (2v/ex) + 4(^2 cos {^^j^ 

^ 1 _ e"2^/v^ 1 + 5^ cos (2v^x) + h cos (4V^x) ' 

and 

e'' =e'^ , (28) 
= - 4e2^/^e^ . (29) 

Thus, we find the relation between the slow-roll parameters for Lp as éP ^ 3(r/'^'^) /4. For 
^2 = 0, we note that is written in terms of ry^^ as 



l-e-2^/v^\[l + 5icos(2Vex)] ^ ^ V^^^i. 

Moreover, rf-'^ has an upper bound as 

&l^4e|5i|. (31) 

In this case, when Ir^^^l reaches the maximum value l^^^axl' minimized. For 62 7^ 0, when 
is minimized, \rf-'^\ does not have to reach |?72iaxl t>ut it is not very different from it. 
Slow-roll infiation occurs for e ^ 1. This corresponds to e"^*^/^ ^ 1 from ( 124|) and ^5^, 
<^ 1 from (I25Ì) . From the slow-roll condition for e*^ we find \!f"P\ <^ 1 automatically. On the 
other hand, the slow-roll condition on "ìf^^ requires \b'2\ <^ 1 for ^ ^ 1, so the x fi^ld gives a 
small contribution to the vacuum energy and infiation is driven dominantly by the radiai mode 
y?. Then, we can also satisfy \é^^\ <^ 1. We note that the slow-roll parameters for x can be 
small in the entire field space. Henceforth, for simplicity, we focus on the case 61 = 6 and 
62 = 0. A non- zero but small 62 would not change the results very much. 
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2.2 Number of e-folds 



Slow-roll inflation ends when e = 1. Since the slow-roll parameter for x remains small, the 
slow-roll condition is violated raainly by ip when q-'^'p^/^ ^ 0.464, with the subscript e denotes 
the end of slow-roll inflation. Here, we refer the readers to later sections for what we here mean 
by "mainly". This gives y^g ~ 0.940. This shows that we always have e"^"^/^ ^ 0.464, which 
is very convenient to compute the number of e-folds and to determine the other field values 

X* and Xe, where -k denotes the moment when the scale of our interest exits the horizon. 

For a product potential, the number of e-folds is found to be [20] 

V 



N 



e^'y,dx. 



(32) 



{V* - Ve) 



(33) 



Thus, from the first equality we obtain 

4 [ Ve 

For = 60, we need e^f*l^ ^ 80.5. Consequently, we can uniquely determine y?^ ~ 5.37. 
From the second equality in (l32l) the evolution of x during slow-roll infiation is governed by 



dN 



1 + ^cos (2Vex) 
27^5 sin {2V^x) 



dx- 



(34) 



Neglecting e ^'/'/v^ -^yj^ich is subdominant, we can analytically integrate this expression to find 



N 



1 
4^5 



log 



tan (v^Xe) 



tan {^/Ix*) 

Thus, ignoring the logarithm multiplied by 5, we obtain 



+ 5 log 



sin {2^Xe) 



sin (2v^x.) 



tan(^V^Xe) ~ ia.n{yix^ 



Thus, for a given x* which may be arbitrary, we can determine Xe uniquely. 
Also note that from Lp^, and (pe-, 



,2f)e 



-2(/3*/V6 
-2<^e/V6 



0.987^ 
0.536. 



so we have 



2he - 2h ^ -0.611 . 
We will use this number in the subsequent sections. 



(35) 



(36) 



(37) 
(38) 

(39) 



3 Perturbations during slow-roll inflation 

In this section, we calculate the primordial perturbations generated during slow-roll infiation 
using the 6N formalism [16j. Because 6N is conformally invariant |21j, we may use the 6N 
formalism safely to compute the primordial perturbations. Note that there do exist isocurvature 
perturbations, and we must follow the post-infiationary evolution e.g. reheating to see how they 
contribute to AT/T and the matter density field. But it is beyond the scope of the present 
work so we do not consider this topic here. 
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3.1 Linear perturbation 

First, let US consider the power spectrum and the spectral index. We can straightforwardly 
find the first derivatives of N with respect to the initial field values (p^, and as [Ì9\ 



e^^-^*^, (41) 



where s* denotes the sign of the potential derivative, i.e. 

+ 1 .f v,>o 

- 1 if < ^ ' 

This is because, while e is positive definite, the derivative of potential may be either positive 
or negative. Then, we can simply substitute (HOj) and ( HTj) into the 6N formulae 

'P( = (f] 7'^'iV,iV,, (43) 



271 J 



ne - 1 = - 2e - ^ ^1 + V^^^^T^ - ^^^^"T^ J ' 

where the right hand sides are evaluated at -k. Here, 7af, is the field space metric with the indices 
a and h being ip and Xi = dN/dcp"' and the same for Ì4, -Raècd is the Riemann curvature 
tensor of the field space, and a semicolon denotes a covariant derivative with respect to the 
field space metric. 

We can conveniently parametrize these quantities from the fact that e = e*^ + e-^. We can 
introduce a dimensionless angle 9 in such a way that [22] 

cos^ 6 =^ , (45) 

smH=—. (46) 
e 

We can see that 6 parametrizes by which field the slow-roll condition is dominated. Further, 
the difference between be and can be fixed from the background as we have seen in ( l39l) . so 
let US write 

2be -2b, = X. (47) 

Then, we can straightforwardly find 

9 V sin^ 9^ 

- 1 = - 2e^ - 4e~2^ 



cos4 (^2 tan2 + tan^ 6 e) 
cos^ é»,, (^ tan^ 9^ - tan^ ^ef , b r, b\ sin^ 6»,, tan^ 9e 

^ ^2 ^nan2 ^, + tan4^e ^"^^ ^ ^ ^an^ + tan^ 

2 , 2 . 2^ tan2 - tan^ 2 (^^ ^an^ ^^rif ^ + tan^ 9eVm , , 

- cos^ 9, tan^ ge .2, 2/^, 4.^ * + ^ .2, -il 4/ ' 49 

9^ + tan* 9e tan^ 6*^ + tan* 9^ 
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-5 -4 -3 -2 -1 
logioé*. 

Figure 1: Logarithmic plot of the boost factor multiplied to the single field power spectrum 
[i/^/(27r)]^ /(2e,t) in f l48p as a function of 6^ and 6e- Fot a sizeable vaine of 6^, this factor may 
be as large as 10^ or even larger. For a large vaine of /nl, however, we may have a moderate 
enhancement of 0{1) ~ (9(10): see Fignre|2]and Table[TJ 

where 

A = [l + (l - e^) tan^ Oe] . (50) 

As we can see, the single fìeld resnlt of "P^ is mnltiplied by a boost factor that depends on 
the valnes of X, 6^ and 6 e- We show the form of this factor as a fnnction of 6^, and 6 e with 
a fìxed X in Figure [H Since this factor is always greater than 1, with the observed value of 
^ 2.5 X 10^^ we may have a considerably lower value of than the usuai single field inflation 
models. Conversely, this factor would suppress the tensor-to-scalar ratio to an unobservable 
level. In our model, the Hubble scale during inflation is given by ~ A/(12^^). When the 
boost factor in the power spectrum is of order unity, we need a/A ~ 10~^^. Therefore, for ^ > 1, 
the self-coupling A must be very small. This is different from the usuai Higgs inflation. 

3.2 Non-Gaussianity 

By extending the SN formalism to second order, we can compute the non-Gaussianity generated 
during slow-roll inflation on super-horizon scales. This corresponds to so-called the locai type 
non-Gaussianity, and is conveniently parametrized by the non-linear parameter /nl [23], which 
is constrained to be —10 < /nl < 74 at 2a level [3|. By taking another derivatives of (HOj) and 



8 



f HT]) . we can obtain [SU] 



dipi 



2e 




3 4fee-4b* 



— + 



XX 



rj_ 
2eì 



-,46e-2fe* 



— + 



with 



+ 1(-1) if y > (<)0. Here, C is given by 



(51) 
(52) 
(53) 



C = 



ere 



^PfiX 



,ere 




where 



^V^jXX _|_ gX^ 



XnW 



(54) 



(55) 



Then, we can compute using the 5N formalism the scale independent part0 of /nl as 



6 ^ 7°'^7'"^iV,,iV,iVrf 

where the right hand side is evaluated at -k. Using the dimensionless angle 9 introduced in (jl5 
and ( H6l) . we can write /nl as 



6. 



-X 



(^2 tan2 + tan^ 



^3 tan^ 6». 



COS^ é'p 



tan^ e. 

cos^ 6*, 



Vie 



y sin^ 61,, tan^ 6^ + tan^ èie 
+ 2e-^ — ^e. 



+ 2 tan^ e. 



cos2 (^2 tan2 + tan^ 
(^tan^é»^ - tan^ O^f 



- sin^ 6^ tan^ tan^ 6»^ 

cos2 (^2 tan2 0^ + tan^ 



ry^^ cos^ df. + sin^ 6'e 



r^r-ee (4cos^^e + ^sin^^, 



(57) 



One can easily verify that among various functions multiplied by the slow-roll parameters, 
the terms proportional to r/^^ allow a huge enhancement at small Q^, and Qe- There are two 
possibilities for large non-Gaussianit}@: 



1. tan2 Q^^ < tan^ 6e and tan2 6^ < tan2 6^. In this case, however, the second terni of (H9Ì) 
may diverge as sin^ 9^/ sin^ 9 e and we need more care. 

2. ^2^an2^^ > tan^6'e and tan2 6'^ < tan2é'e. In this case, r/'s give dominant contributions 
to n^. Therefore we focus on this case for the discussion below. 



^We do have a scale dependent part of /nl, which comes from the intrinsic non- Gaussiani ty in the field 
fluctuations. But it is always inuch sinaUer than 1 [24^ so we do not consider it here. 
^The two conditions are closely related to those in Ref. [22] . 
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Figure 2: The prefactor of the last term of (1571) as a function of 9^, and Of.. There is a region 
where this factor blows up so that we can obtain a detectable value of |/nl|- Note that we may 
have even larger regions that allow detectable /nl for X > 0, which however means <y9^ < ipf. so 
we do not consider this case. 

Note that for À^tan^Oi, > tan^6'e and tan^ 6*^ > tan^ we have a very small value of /nl, 
consistent with current observations though, and we do not consider this case. For J\? tan^ 9^, < 
tan^ 9e and tan^ 9i, < tan^ 6'e, /nl could be large. But this is not possible because > e*^ always, 
which means inflation occurs mostly along x direction. 

In Figure [21 we show the prefactor of the last term of (1571) as a function of 9^ and 9e with 
X fixed by (!39|) . The condition for large /nl is clear from this plot: the contribution of to 
the slow-roll condition should be larger at the end of inflation than at the moment of horizon 
Crossing to occupy 10""^% ~ 10% of e^. This means along the x direction, the situation is more 
or less similar to the hilltop inflation, where the initial condition of placing x* dose to the locai 
maximum of its cosine potential is not as severe as it sounds [25J. 

We can estimate the value of /nl as foUows. For 6'e -C 1, and for J\?tSiV?9i, > tan^6'e and 
tan^ 9i, < tan^ 6'e, from the ì]^^ terms of fl57j) . we fìnd 



b f _x tan'^6'e vv , ^ • 2 /, tan'^6'e \ 



(58) 



From fl30|) . in the region with a negative i]'^^^ using e 2v*/v^ ^ ]^ ^-^^ g 1h>^I^ ^ \ 
replace r/^'^'s by 



/2, we can 




(59) 



(60) 
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10 0.00340 10"^ 10-2-1^ 0.952 -14.6 -16.5 
1 0.0380 10-5 10-2-^5 0.959 -60.2 -63.7 



Table 1: and /nl with given 6*^, 6'e, ^ and S. The vaine of "P^ is set by the COBE normahzation 
2.5 X 10^^. Note that for these valnes of 9-^, 9e and ^, we find the boost factor of the power 
spectrum 1.06 and 1.03 respectively. In the last column we show the analytic estimate of /nl 
given by fl58l) . 

It is also possible that r]}^ is positive for C,\6\ > 0.1 from fl36|l for x* being dose to a locai 
maximum. But, in this case, rj^^ ~ —0.1, so the spectral index would be modifìed too much by 
X from (jin]). Therefore, for a typical vaine of \ri^^\ ~ 0.01 during inflation, we restrict ourselves 
to the region with r]}^ < 0. Note that there are inequalities between 9^,, 9 e and ^S"^ as 



Thus, by choosing appropriate valnes of 9^ and ^ in (!58|) and taking into account the evolution of 
the X field from (!36l) for the corresponding vaine of 9e , we can have large /nl with observationally 
consistent n^. In Table [1] we show two sets of sample valnes. Note that /nl < in models 
where non-Gaussianity is generated by a divergent trajectory falling of a ridge [221 126] . 



To summarize, we have studied a simple multi-field inflation model by introducing a complex 
scalar fìeld with non-minimal coupling to gravity. We have considered the general action for 
the complex scalar fìeld with dimension-4 interactions in the potential. Assuming that the 
non-minimal coupling breaks the global U{1) symmetry by a small amount, the Lagrangian has 
been brought to a simple form with a non-canonical kinetic term for the angular mode x ci^nd 
a product form of the potential for x and the modulus ip. Then, the potential of ip becomes 
fiat at large fìeld valnes due to the non-minimal coupling as in Higgs inflation, while a small 
breaking of the global U{1) symmetry in the potential terms makes the potential for x Aa-t 
as well. This allows slow-roll two-fìeld inflation. Using the 6N formalism, we have computed 
the power spectrum P^, the spectral index ra^ and the non-linear parameter /nl- With the 
hilltop type initial condition of x, we have shown that /nl can be large enough to be detected 
in near future. The generality of the hilltop type inflation suggests that the initial conditions 
for large non-Gaussianity are not as fìnely tuned as previously believed, and the possibility of 
large enough /nl is rather high. We have also presented sample sets of the parameters which 
give |/nl| = C(10) as well as the observationally consistent n^. 

Our toy model with a complex scalar fìeld can be easily embedded in well-known particle 
physics models such as the Standard Model with two Higgs doublets. When the global U{1) 
symmetry adopted in this paper is identifìed with the Peccei-Quinn (PQ) symmetry, in the 
limit of a small violation of the PQ symmetry the pseudo-scalar Higgs boson can obtain a fiat 




(61) 



sin^ 9e <4:^5^ . 



(62) 



4 Conclusions 
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potential such that it is an inflaton candidate. It will be interesting to see what the constraints 
on two-field inflation wc havc found in this paper imply for the scarch of the additional Higgs 
boson at the Large Hadron Colhder. We leave this work in a future study. 

Acknowledgements 

We thank Ki Young Choi for helpful comments. JG is grateful to the Institute for Advanced 
Study at the Hong Kong University of Science and Technology for hospitahty while this work 
was under progress. This work is partiaUy supported by Korean-CERN feUowship. 



12 



References 



A. H. Guth, Phys. Rev. D 23, 347 (1981) ; A. D. Linde, Phys. Lett. B 108, 389 (1982) ; 
A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982). 

See e.g. A. R. Liddle and D. H. Lyth, "Cosmologica! inflation and large-scale structure," 
Cambridge, UK: Univ. Pr. (2000) 400 p ; V. Mukhanov, "Physical foundations of cosmol- 
ogy," Cambridge, UK: Univ. Pr. (2005) ^21 p ; S. Weinberg, "Cosmology," Oxford, UK: 
Oxford Umv. Pr. (2008) 593 p. 

E. Komatsu et al. [WMAP Collaboration] , Astrophys. J. Suppl. 192, 18 (2011) 
[a rXiv:1001.4538 [astro-ph.CO]]. 

For a review, see e.g. D. H. Lyth and A. Riotto, Phys. Rept. 314, 1 (1999) 
|a rXiv:hep-ph/9807278j . 



[Planck Collaboration], arXiv:astro-ph/0604069, 



N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 103 (2004) 
|arXiv:astro-ph / 0406398 1 ; Class. Quant. Grav. 27, "Focus section on non-linear and non- 
Gaussian cosmologica! perturbations" (2010) ; Adv. Astron. 2010, "Testing the Gaussian- 
ity and Statistica! Isotropy of the Universe" (2010). 

T. Tanaka, T. Suyama and S. Yokoyama, Class. Quant. Grav. 27, 124003 (2010) 
^arXiv:1003.5 057 [astro-ph.CO]]. 

B. L. Spokoiny, Phys. Lett. B 147, 39 (1984) ; F. S. Accetta, D. J. Zoller and M. S. Turner, 
Phys. Rev. D 31, 3046 (1985) ; D. S. Salopek, J. R. Bond and J. M. Bardeen, Phys. Rev. 
D 40, 1753 (1989). 

F. L. Bezrukov and M. Shaposhnikov, Phys. Lett. B 659, 703 (2008) (arXiv:0710.3755l 
[hep-th]]. 

R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH Collaboration 
and DELPHI Collaboration and L3 Collaboration and OPAL Collaboration], Phys. Lett. 
B 565 (2003) 61 [ar Xiv:hep-ex/ 0306033|. 

C. P. Burgess, H. M. Lee and M. Trott, JEEP 1007 (2010) 007 [arXiv: 1002.27301 [hep-ph]] . 

G. F. Giudice and H. M. Lee, Phys. Lett. B 694, 294 (2011) [arXiv:1010.1417l [hep-ph]]. 

C. P. Burgess, H. M. Lee and M. Trott, JHEP 0909 (2009) 103 [iiXw: 0902. 4465 [hep- 
ph]]; J. L. F. Barbon and J. R. Espinosa, Phys. Rev. D 79 (2009) 081302 (2009) 081302 
[arXiv:0903.0355 [hep-ph]]: M. P. Hertzberg, a rXiv:1002.2995 [hep-ph]. 

F. Bezrukov, A. Magnin, M. Shaposhnikov and S. Sibiryakov, JHEP 1101 (2011) 016 
|arXiv:1008.5T57 ' [hep-ph]]. 

R. N. Lerner and J. McDonald, Phys. Rev. D 82 (2010) 103525 [arXiv: 1005.2"978l [hep-ph]] . 



13 



[16] A. A. Starobinsky, JETP Lett. 42 (1985) 152 [Pisma Zh. Eksp. Teor. Fiz. 42 
(1985 ) 124] ; M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996) 
|arXiv:astro-ph/9507001| ; M. Sasaki and T. Tanaka, Prog. Theor. Phys. 99, 763 (1998) 
|arXiv:gr-qc/98010l7| ; J. O. Gong and E. D. Stewart, Phys. Lett. B 538, 213 (2002) 
|arXiv:astro-ph/0202098| ; D. H. Lyth and Y. Rodriguez, Phys. Rev. Lett. 95, 121302 
(2005) |arXiv:astro-ph/0504045|. 

[17] H. M. Lee and A. Papazoglou, Phys. Rev. D 80 (2009) 043506 [a rXiv: 0901. 49621 [hep-th]]. 

[18] H. Y. Chen and J. O. Gong, Phys. Rev. D 80, 063507 (2009) [arXiv:0812.4649 [hep- 
th]] ; H. Y. Chen, J. O. Gong, K. Koyama and G. Tasinato, JCAP 1011, 034 (2010) 
[a rXiv: 1007.2068 [hep-th]]. 

[19] J. Garcia-Belhdo and D. Wands, Phys. Rev. D 53, 5437 (1996) [arXiv: astro-ph/ 95 11029]. 

[20] K. Y. Choi, L. M. H. Hall and C. van de Bruck, JCAP 0702, 029 (2007) 
(arXiv:astro-ph/Ó701247| . 

[21] N. Makino and M. S asaki, Prog. Theo r. Phys. 86, 103 (1991) ; S. Koh, J. Korean Phys. 
Soc. 49, S787 (2006) |arXiv:astro-^h70510030| ; T. Chiba and M. Yamaguchi, JCAP 0810, 
021 (2008) iarXiv: 0807^ 965 [astro-ph]] ; N. Sugiyama and T. Futamase, Phys. Rev. D 81, 
023504 (2010) ; J. O. Gong, J. c. Hwang, W. I. Park, M. Sasaki and Y. S. Song, JCAP 
1109, 023 (2011) [ar Xiv: 1107. 1840 [gr-qc]]. 

[22] C. T. Byrnes, K. Y. Choi and L. M. H. Hall, JCAP 0810, 008 (2008) ;arXiv:0807.lT0T] 

[astro-ph]]. 



[23] E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 (2001) [ arXiv:astro-ph/0005036| . 
[24] F. Vernizzi and D. Wands, JCAP 0605, 019 (2006) [arXiv:astro-ph/0603799| . 



[25] L. Boubekeur and D. H. Lyth, JCAP 0507, 010 (2005) [ar Xiv:hep -ph/0502047] ; K. Kohri, 
C. M. Lin and D. H. Lyth, JCAP 0712, 004 (2007) [a rXiv:0707.3826 [hep-ph]]. 

[26] J. EUiston, D. J. Mulryne, D. Seery and R. Tavakol. larXiv:1106.2153l [astro-ph. CO]. 



14 



